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1.  Foreword 

The  following  document  provides  a  final  progress  report  for  the  work  completed  under  a  Short 
Term  Innovative  Research  (STIR)  award  from  the  US  Army  Research  Office  to  a  team  from  The 
Pennsylvania  State  University  (PSU)  and  Georgia  Institute  of  Technology  (GT).  The  Pis  for  this 
project  were  C.  Saldana  (GT)  and  R.  Voigt  (PSU).  The  STIR  project  was  initiated  at  PSU  with  C. 
Saldana  as  PI,  whom  subsequently  transferred  primary  employment  to  GT  during  the  contracting 
period.  The  project  was  then  subcontracted  to  GT  by  PSU. 
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4.  Statement  of  the  Problem  Studied 

The  focus  of  the  present  study  was  on  establishing  principles  to  support  multi-scale  surface 
texturing  of  metallic  components  for  use  in  military  performance  applications.  In  this  regard,  two 
different  surface  texturing/modification  methods  were  considered:  Surface  Severe  Plastic 
Deformation  (S2PD)  and  Modulation-Assisted  Texturing  (MAT).  For  S2PD,  the  effects  of 
processing  parameters  on  the  resulting  subsurface  deformation  field  in  terms  of  both  mechanics 
and  micro  structure  (e.g.,  micro  structure,  texture)  were  quantified  for  a  unit  model  process  as  well 
as  for  an  aggregate-level,  strain-induced  surface  modification  configuration.  For  MAT,  the  effects 
of  surface  texture  form  (e.g.,  micro-scale  topography)  on  surface  integrity  measures  and 
tribological  wear  performance  were  quantified.  The  ensuing  results  are  useful  as  a  preliminary 
step  toward  realizing  a  hybrid  model-based  processing  framework  for  designing  multi-functional 
surfaces  with  significantly  enhanced  mission-specific  performance  capabilities. 
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5.  Summary  of  the  Most  Important  Results 

The  primary  project  results  are  summarized  below  in  terms  of  experimental  and  modeling  results 
associated  with  the  surface  severe  plastic  deformation  and  controlled  surface  texturing  methods 
investigated. 

5.1  Surface  Severe  Plastic  Deformation  (S2PD) 

In  the  present  project,  the  evolution  of  microstructure  and  crystallographic  texture  during  S2PD  in 
SMAT  were  investigated.  As  the  microstructure  created  during  SMAT  is  a  consequence  of  the 
coupled  effect  of  multiple  surface  impacts/indentations,  insights  were  assembled  from  elucidating 
the  mechanics  and  microstructure  evolution  in  a  unit  surface  indentation  process  (e.g.,  unit  model 
mechanics)  and  using  these  models  to  establish  effects  in  a  scaled  configuration  (e.g.,  S2PD 
processing). 

5.1 .1  Unit  Model  Mechanics 

Digital  Image  Correlation  (DIC)  was  used  to  characterize  material  flow  during  surface 
indentation.  For  this,  the  face  of  the  workpiece  (plane  Y-Z  in  Fig.  1)  exposed  through  the 
tempered  glass  window  was  illuminated  during  indentation  using  a  bright  light  source,  which 
resulted  in  a  concentrated  speckle  pattern  on  its  surface.  Material  flow  was  recorded  in  a  sequence 
of  digital  images  captured  using  a  high-speed  camera.  Subsequently,  movement  of  these  speckles 
due  to  deformation  imposed  during  indentation  was  delineated  using  image  correlation 
algorithms.  Subroutines  were  thereafter  implemented  to  extract  several  deformation  parameters 
including  strain  rate  fields,  strain  fields  and  rotation  fields.  To  do  this,  a  grid  was  defined  in  the 
region  of  interest  under  the  indenter.  Subsequently,  path  lines  followed  by  each  point  in  the  grid 
were  extracted  and  deformation  experienced  by  the  corresponding  point  in  the  grid  were 
integrated  along  the  path  line.  In  this  manner,  deformation  histories  of  different  material  points 
under  the  impacting  media  were  characterized.  The  computational  routines  were  written  in-house 
using  Matlab  [1].  The  deformation  histories  were  coupled  with  the  Visco-Plastic  Self-Consistent 
polycrystal  (VPSC)  framework  to  simulate  evolution  of  crystallographic  textures  [2] . 


Figure  1.  In-situ  characterization  of  material  flow  during  surface  deformation.  Characterization  of  (a) 
deformation  fields  to  produce:  (b)  effective  strain,  (c)  rotation  and  (d)-(f)  Lagrangian  strain  components. 

Orientation  imaging  microscopy  (OIM)  of  deformed  samples  was  performed  using  electron 
backscatter  diffraction  (EBSD).  Specimens  were  prepared  by  mounting  deformed  samples  in 
epoxy,  followed  by  mechanical  polishing  to  0.04  um  finishes.  Final  mechanical  polishing  was 
performed  using  colloidal  silica  for  longer  periods  for  enhancing  indexing.  Subsequently,  samples 
were  also  polished  with  a  Hitachi  4000  plus  ion  beam  polisher  for  5  mins.  Scan  step  sizes 
between  0.1  um-1  um  were  used.  Larger  scan  step  sizes  facilitated  quick  collection  of  statistically 
significant  data  for  studying  crystallographic  textures.  Beam  diameter  was  set  at  150  nm  and 
accelerating  voltage  was  set  at  20  KV.  Analysis  of  crystallographic  textures  extracted  from  EBSD 
micro-graphs  was  performed  using  MTEX  software.  Inverse  pole  figures  (IPFs)  illustrating 
empirically  observed  evolution  of  microstructure  at  different  locations  under  the  indenter  are 
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shown  in  Fig.  2.  The  grains  exhibit  a  pancaked  shape  closer  to  the  indenter-workpiece  interface 
due  to  shear  imposed  (region  I).  However,  this  effect  decays  with  strain  away  from  the  indenter, 
as  evidenced  from  near  uniaxial  grain  morphologies  observed  in  these  regions  (Fig.  2,  regions  II 
and  III). 


Figure  2.  Orientation  imaging  microscopy  underneath  unit  surface  deformation  test.  Three  regions  (I,  II  and 
III)  were  chosen  for  performing  OIM.  Refer  inverse  pole  figure  inset  for  color  code. 

Crystallographic  textures  comprising  orientation  distribution  functions  (ODFs)  and  pole  figures 
(Fig.  3)  were  extracted  from  these  IPFs  by  discrete  binning  of  the  orientations.  Significant 
differences  between  crystallographic  textures  extracted  from  the  three  different  locations  were 
seen,  owing  to  the  heterogeneous  deformation  under  the  indenter  manifesting  through  different 
strain  tensor  components  (Fig.  1).  The  strain  tensor  components  were  calculated  using  F  = 

-(Fr  ■  F  —  /)  where  I  is  the  identity  tensor.  For  instance,  crystallographic  textures  in  region  I 
(Fig.  3)  exhibit  a  dominant  simple  shear  character  featuring  concentrations  along  fibers 
reminiscent  of  the  fl5  f2  and  f3  fibers  conventionally  seen  in  simple  shear  deformation  geometries 
[3].  The  Lagrangian  strain  tensor  components  in  this  region  were  En  =  0.1,  E22  =  -0.1  and  E12  = 
0.9.  En  and  E22  refer  to  the  horizontal  EYY  and  vertical  Ezz  Lagrangian  strain  components, 
respectively.  On  the  other  hand,  crystallographic  textures  in,  region  II  (Fig.  3)  exhibit  a  uniaxial 
tension  character,  with  the  tensile  axis  oriented  along  the  horizontal  (Y)  direction  [4].  This  is  due 
to  the  dominant  uniaxial  strain  component  in  region  II  with  En  =  0.2,  E22  =  -0.2  and  E12  =  0.01 . 

Evolution  of  crystallographic  textures  during  indentation  was  modeled  using  the  VPSC 
framework.  For  this,  the  deformation  histories  were  obtained  from  in  situ  DIC  based 
characterization  of  the  deformation  zone  during  indentation.  Pre-deformation  crystallographic 
textures  were  characterized  by  OIM  using  EBSD  and  revealed  lack  of  preferred  orientations  (e.g., 
random  crystallographic  textures),  as  in  Fig.  3.  Hence,  the  pre-deformation  discretized  ODFs  for 
input  to  the  VPSC  framework  were  produced  by  sampling  from  uniformly  distributed  orientation 
space.  Numerically  simulated  crystallographic  textures  in  surface  indentation  are  shown  in  Fig.  3. 
Here,  parameters  of  the  VPSC  framework  were  tuned  until  a  good  match  between  simulated  and 
empirically  obtained  results  was  found.  It  was  seen  that  a  grain  fragmentation  scheme  with  aspect 
ratio  R=5  was  adequate  for  simulating  the  texture  evolution.  This  scheme  of  grain  fragmentation 
is  in  the  same  vein  as  mechanics  governing  geometric  dynamic  recrystallization,  which  involves  a 
single  long,  pancake- shaped  grain  pinching  off  into  several  small  grains  upon  meeting 
mechanistically  defined  aspect  ratios  [5].  While  some  discrepancies  were  seen  in  the  maximum 
intensities  of  the  empirical  and  simulated  crystallographic  textures  (Fig.  3),  the  calibrated  VPSC 
frameworks  essentially  reproduced  all  qualitative  aspects  of  the  post-deformed  crystallographic 
textures.  Variations  in  intensities  arose  due  to  large  grain  sizes  manifesting  in  concentrations  in 
the  orientation  space. 
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Figure  3.  Empirical  and  simulated  crystallographic  textures  corresponding  to  subsurface  regions  after 
surface  indentation.  (Ill)  and  (022)  pole  figures  as  well  as  02  sections  (0°  and  45°)  are  shown. 

5.1 2  S2PD  Processing 

Surface  mechanical  attrition  treatment  (SMAT)  was  performed  on  as-received  OFHC  Cu  (alloy 
101:  MSC  direct,  2.3  mm  thick).  For  SMAT,  spherical  silica  Si02  shots  (diameter  =  200  um) 
were  accelerated  using  a  directed  blast  of  air  at  90  Psi  (0.4  MPa).  In  this  manner,  the  surface  was 
treated  with  several  hundred  impacts  oriented  perpendicularly  with  respect  to  the  surface  normal 
(Z  axis  in  Fig.  1)  whereby  directionality  of  SMAT  was  maintained.  SMAT  was  performed  for 
600  s  and  1200  s.  The  velocity  of  impact  was  measured  to  be  approximately  10  m/s  by  recording 
the  mass  flow  rate.  Numerical  simulation  was  performed  using  ABAQUS  6.14  as  a  sequence  of 
impacts,  implemented  through  in-house  python  subroutines  using  the  scripting  interface  in 
ABAQUS.  The  mechanical  state  of  the  workpiece  prior  to  impacts  were  extracted  from  the  post¬ 
deformation  state  left  after  the  previous  impact  simulation.  Python  subroutines  were  used  to 
achieve  this  feature  whereby  the  deformation  history  at  the  end  of  an  impact  is  not  lost  and 
material  response  analogous  to  that  in  real  SMAT  comprising  a  sequence  of  impacts  is  simulated. 
Impacting  shots  were  modeled  as  spherical  elastic  bodies  with  (d  =  200  um,  E  =  300  GPa,  v  = 
0.21).  The  workpiece  was  modeled  as  an  isothermal  body  (C3D84  elements)  with  reduced 
integration  and  hourglass  control.  The  iso-thermal  state  of  the  workpiece  was  governed  by  the 
diffuse  nature  of  deformation  imposed  by  spatially  and  temporally  discrete  impacts  across  the 
surface  of  the  workpiece,  which  prohibit  concentrated  heat  dissipation  and  temperature  rise  due  to 
plastic  deformation.  The  workpiece  was  modeled  as  a  Johnson-Cook  material:  o  —  (A  + 

(rr _ rr  \  TTL 

1 - y-j  ,  where  [A,  B,  C,  £0,  n,  m]  =  [90  MPa,  292  MPa,  0.025,  1  s'1, 

T-m  TrJ 

0.31,  1.09],  respectively.  Gaussian  distributed  velocities  were  assigned  to  the  impacting  shots. 
Post  simulation  processing  was  performed  to  couple  data  obtained  from  numerical  FEM  based 
simulations  with  the  VPSC  framework  for  simulating  evolution  of  crystallographic  textures 
explicitly.  This  was  done  using  Python  subroutines,  whereby  displacement  undergone  by  element 
nodes  close  to  the  surface  was  converted  to  velocity  gradient  tensors,  which  comprises  the  raw 
input  for  the  VPSC  framework. 

Figure  4  shows  the  response  of  material  to  an  individual  impact  V  =  -0.5ex+0.6eY-10ezm/s.  Here 
ex,  eY  and  ez  are  unit  vectors.  Lagrangian  strain  components  reveal  a  3D  material  flow  pattern 
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involving  comparable  strain  components  featuring  max  IEnl  within  the  same  order  of  magnitude. 
Here  IEnl  refers  to  components  in  the  Lagrangian  strain  tensor  E_.  Implications  of  the  3D  nature  of 
the  deformation  field  lie  in  strain  path  changes,  reversals  and  spatial  strain  heterogeneity  in  the 
final  strain  fields.  These  have  further  structural  and  functional  manifestations  in  the  resulting 
product,  such  as  accelerated  micro  structure  refinement  during  SPD  [6,  7].  A  maximum  effective 
strain  rate  edot  =  5.3  x  105  s'1  was  seen  under  the  impacting  shot.  Figure  5  shows  the  numerically 
simulated  strain  profiles  under  the  surface  treated  using  SMAT.  These  curves  elucidate 
representative  mean  strain  profiles  obtained  from  several  different  locations  under  the  surface. 
Interestingly,  a  constant  width  of  material  accommodates  the  strain  imposed  during  SMAT, 
irrespective  of  the  number  of  impacts  on  the  surface.  Despite  a  complicated  deformation 
geometry  ensuing  from  the  coupled  effect  of  several  impacts,  the  constancy  in  dimensions  of  the 
affected  zone  is  reminiscent  similar  observations  for  deformation  in  burnishing  and  machining. 
The  simulation  was  terminated  at  the  end  of  850  impacts,  at  which  point  a  peak  effective  strain  8 
=  1.3  was  seen  numerically.  The  strain  field  exhibited  a  near-Gaussian  spatial  distribution  at  the 
end  of  the  simulation,  arising  from  Gaussian  distributed  velocities  of  the  impacting  shots. 
Additionally,  this  behavior  was  coupled  with  an  exponential  decay  in  magnitude  under  the 
surface,  similar  to  machining,  despite  different  underlying  deformation  geometries. 
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CONFIGURATION  LAGRANGIAN  STRAIN  COMPONENTS  (x  10  2) 

Figure  4.  Simulated  response  of  material  to  a  single  surface  impact.  Configuration  of  the  numerical  model 
is  shown  on  the  left.  Effective  and  Lagrangian  strain  tensor  field  components  are  shown  on  the  right. 
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Figure  5.  Numerically  obtained  strain  profiles  under  the  surface  deformed  by  SMAT.  The  inset  shows  the 
strain  field  after  850  impacts  and  a  cross  sectional  view  illustrating  the  strain  profile  under  the  surface. 


Characterization  of  crystallographic  textures  in  as-received  and  post-SMAT  samples  was 
performed  in  a  Panalytical  X'Pert  pro  MRD  x-ray  diffractometer  equipped  with  a  Cu  x-ray 
source,  an  Eulerian  cradle  goniometer  and  a  solid  state  X'celerator  detector.  Measurements  were 
performed  at  45  KV  and  40  mA  using  the  Kal  radiation  of  Cu  k  =  1 .540598  angstroms  with  a  scan 
step  size  of  5  degrees.  For  this,  the  surface  perpendicular  to  the  Z  direction  was  exposed  to  the 
incident  X  ray  beam.  Pole  figures  corresponding  to  the  first  three  peaks  in  Cu  (111),  (200)  and 
(220)  at  diffraction  angles  of  43,  50  and  74  degrees,  were  measured  respectively.  The  diffraction 
angles  were  fine  tuned  before  making  each  measurement.  Post  measurement  analyses  including 
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pole  figure  inversion  of  crystallographic  textures  were  performed  using  MTEX.  Figure  6  shows 
ODFs  and  (111)  and  (022)  pole  figures  of  crystallographic  textures  obtained  from  surfaces,  pre- 
SMAT  and  post-SMAT  treatments  (600  s  and  1200  s).  Prior  to  SMAT,  the  surface  exhibits  a  cube 
texture  with  pole  figures  featuring  2/m  symmetry.  Upon  SMAT,  the  pole  figures  assumes  360- 
degree  rotation  symmetry  about  the  central  axis  that  corresponds  to  normal  direction  with  respect 
to  the  surface.  Additionally,  crystallographic  textures  obtained  from  600  s  and  1200  s  samples 
show  negligible  variation  with  respect  to  each  other  suggesting  a  saturated  state  of  the  underlying 
microstructure  with  respect  to  the  deformation  imposed  during  SMAT.  Such  saturated  micro 
structural  states  are  commonly  seen  during  imposition  of  SPD  [8]. 


(Ill)  (022)  (111)  (022) 


Figure  6.  Empirical  and  simulated  (111),  (022)  pole  figures  from  initial  and  post-SMAT  samples. 

Numerically  obtained  material  flow  during  SMAT  was  integrated  with  the  VPSC  framework 
explicitly  to  simulate  evolution  of  crystallographic  textures.  This  was  done  by  extracting  the 
velocity  gradient  tensors  from  arbitrarily  chosen  elements  directly  under  the  surface  of  the 
numerically  simulated  workpiece.  To  this  end,  elements  close  to  the  center  of  the  workpiece 
exhibiting  peak  in  the  Gaussian  distributed  strain  field  were  chosen.  After  subsequent  post¬ 
processing,  a  compiled  deformation  history  file  containing  material  response  data  from  the 
impacts  was  produced  and  input  to  the  VPSC  framework.  Here,  deformation  steps  resulting  in  an 
incremental  effective  strain  were  ignored  assuming  a  purely  elastic  response.  Such  events  are 
manifestations  of  elastic  waves  from  impacts  in  non-neighboring  regions  of  the  arbitrarily  chosen 
element.  The  starting  un-deformed  texture  which  also  comprises  input  to  the  VPSC  framework  in 
a  discretized  ODF  form  (i.e.,  list  of  grain  orientations  corresponding  to  the  ODF)  was  obtained  by 
statistically  sampling  grains  from  the  empirically  measured  pre-deformed  ODF  (Fig.  6). 

Simulated  crystallographic  textures  (ODFs,  pole  figures)  are  shown  in  Figs.  6  and  7,  respectively, 
and  exhibit  qualitative  similarity  with  the  empirical  counterparts.  This  is  evidenced  from  the 
numerically  simulated  ODFs  which  exhibit  concentrations  along  fibers  featuring  02  symmetry 
(black  arrows  in  Fig.  7),  which  correspond  to  the  360°  symmetry  of  the  deformed  sample.  The 
location  of  these  fibers  was  coincident  with  those  obtained  empirically  through  600  s  and  1200  s 
of  SMAT.  The  discontinuous  nature  of  the  fibers  obtained  numerically  is  a  consequence  of 
spatially  discontinuous  deformation  imposed  during  the  850  impacts.  However,  this  was  not  the 
case  empirically  owing  to  the  potentially  larger  number  of  impacts  that  result  in  a  more 
continuous  fibers  distributed  through  the  orientation  space.  In  order  to  understand  these  effects, 
360°  02  symmetry  was  imposed  on  the  numerically  obtained  crystallographic  textures, 
representative  of  the  numerically  obtained  saturated  state  of  the  workpiece  (Fig.  4).  Results 
obtained  using  this  approach  coincide  well  with  the  empirically  obtained  crystallographic 
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textures,  as  seen  in  comparison  of  the  pole  figures  as  well  as  ODFs.  On  further  analysis,  it  was 
seen  that  the  crystallographic  textures  obtained  during  SMAT  share  similarities  with  compression 
textures,  which  constitute  fibers  coincident  with  those  obtained  empirically  during  SMAT. 
Reasons  for  this  similarity  constitute  strain  path  changes  involving  reversals  that  negate  other 
Lagrangian  strain  components  that  resulted  from  previous  impacts  and  additively  amplify  the 
compression  E33  component  under  the  surface  during  SMAT.  This  elucidation  has  consequences 
for  realizing  process-structure-function  frameworks  through  micro  structure  engineering  using 
SMAT.  The  numerical  framework  essentially  enables  material  flow  field  prediction  as  a  function 
of  parameters  of  the  SMAT  process,  which  include  impacting  shot  geometry,  shot  density, 
velocity  and  angle  of  incidence. 


Figure  7.  Empirical  (600  s,  1200  s)  and  simulated  (850  impacts,  360  02  symmetry)  orientation  distribution 
functions.  Black  arrows  show  coincidence  of  dominant  fibers  in  empirically  obtained  and  numerically 
simulated  samples.  02  sections  =  0°  and  45°  are  shown. 


Results  from  these  experimental  and  calibrated  numerical  investigations  can  be  directly  integrated 
within  the  expected  functional  response  space  in  order  to  optimize  surface  engineering  processes 
involving  SMAT.  The  framework  described  is  directly  relevant  to  a  range  of  bio-medical  implant 
engineering  technologies,  especially  those  involving  anisotropically  brittle  materials,  such  as  Mg 
alloys.  For  instance,  anisotropy  inhibits  formability  and  machinability  and  prohibits  widespread 
use  of  Mg  alloys  in  bio-medical  applications.  As  in  the  above,  surface  engineering  technologies 
such  as  SMAT  can  be  used  to  tune  the  crystallographic  textures  of  Mg  using  this  framework  in 
order  to  engineer  starting  microstructures  for  engendering  higher  formability  and  machinability. 


5.2  Multi-Scale  Texturing 

In  the  present  project,  micro-scale  textures  of  controlled  size  and  morphology  were  generated  by 
control  of  underlying  process  parameters  in  modulation-assisted  texturing.  Modulation-assisted 
texturing  (MAT)  broadly  refers  to  a  class  of  surface  generation  processes  that  utilize  controlled 
modulation  of  the  workpiece  or  tool  during  conventional  machining.  Two  modulation 
configurations  are  considered  wherein  the  superimposed  modulation  is  applied  in  the  direction  of 
the  workpiece  centerline.  The  first  configuration  is  that  of  a  facing  process  with  modulation 
applied  in  the  depth  of  cut  direction,  here  referred  to  as  plunging-type  texturing  and  depicted  in 
Fig.  8(a).  In  this  configuration,  the  tool  is  engaged  at  a  constant  depth  of  cut  (d)  along  the 
workpiece  axis  and  traverses  radially  along  the  end  face  at  a  constant  feed  rate  (s0).  The 
modulation  in  this  case  is  applied  in  the  depth  of  cut  direction,  which  causes  a  periodic  change  in 
depth  of  cut.  Depending  on  the  selected  parameters,  the  tool  is  periodically  separated  from  the 
workpiece,  causing  discretization  of  the  chip  formation  and  yielding  a  periodic  texture  on  the 
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workpiece  face.  A  second  configuration  is  that  of  cylindrical  turning  with  modulation  applied  in 
the  feed  direction,  here  referred  to  as  sliding-type  texturing  and  depicted  in  Fig.  8(b).  In  this 
configuration,  the  tool  is  engaged  at  a  constant  depth  of  cut  (d)  and  travels  parallel  to  the  axis 
along  the  part  length  at  a  constant  feed  rate  (s0).  The  modulation  is  applied  in  the  direction  of  the 
undeformed  chip  thickness  and  the  interaction  of  successive  machining  passes  imparts  a  periodic 
micro-scale  topography  to  the  surface.  Micro-scale  depression-type  textures  (dimples)  were 
generated  using  these  MAT  processes  and  corresponding  surface  models  were  developed  as  a 
framework  for  their  manufacture.  The  tribological  performance  of  the  dimple-type  textures  was 
characterized  by  ball-on-flat  wear  testing.  The  effect  of  the  texturing  process  on  surface  and 
subsurface  microstructure  was  measured  using  EBSD. 


Figure  8.  Schematic  showing  tooling  configuration  to  produce  textured  surfaces  in  (a)  plunging-type 
texturing  (facing)  and  (b)  sliding-type  texturing  (turning). 


5.2.7  Surface  Modeling 

A  geometric  surface  model  for  plunging-type  texturing  can  be  determined  by  considering  the 
interaction  of  the  tool  with  the  work  surface.  For  facing  configurations  in  the  absence  of 
superimposed  modulation,  motion  of  the  texturing  tool  is  restricted  to  a  plane  parallel  to  the  work 
surface  and  takes  the  shape  of  an  Archimedan  spiral  (i.e.,  r  =  a  +  0b,  a  and  b  are  constants)  as  in 
Eq.  (1), 


2*71 


Eq.  (1) 


where  s0  is  the  feed  rate,  r  the  radial  coordinate  and  0  the  angular  coordinate  as  in  Fig.  9(a). 
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Figure  9.  Tool  motions  in  (a)  plunging-type  texturing  (facing)  and  (b)  sliding-type  texturing  (turning). 
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In  the  case  of  constant  angular  frequency  (i.e.,  constant  fw),  d0/dt  is  constant  and  the  relative 
phase  angle  between  successive  machining  passes  [9,  10]  is  given  by  Eq.  (2), 


cp  = 


Eq.  (2) 


In  Eq.  (2),  in-phase  modulation  conditions  correspond  to  cp  =  n*27i  and  out-of-phase  modulation 
conditions  correspond  to  cp  =  (2n+l)7i,  where  n  is  zero  or  a  non-zero  integer.  Intermediate  phase 
differences  result  for  other  values  of  cp.  For  constant  angular  frequency,  a  relationship  between 
tool  position  (r,  0,  z)  and  time  t  is  relatively  straightforward.  Constant  surface  speed 
configurations  are  important  for  facing-type  machining  configurations  as  such  conditions  provide 
for  uniform  material  removal  conditions  (e.g.,  cutting  velocity)  and  also  yield  surface  textures  of 
constant  size,  as  is  shown  in  the  ensuing.  In  this  case,  d0/dt  is  time-varying  to  maintain  constant 
v0  and  the  distance  traveled  along  the  spiral  must  be  considered.  The  path  length  for  a  given 
angular  traverse  (i.e.,  to  position  0i)  is  given  by  the  arc  length  of  the  spiral  in  Eq.  (3), 
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Eq.  (3) 


where  02  =  27iR/s0  is  the  angular  position  (i.e.,  number  of  spiral  turns)  at  the  periphery  of  the 
workpiece  face  and  R  is  the  radius  of  the  workpiece  (Fig.  9).  At  constant  surface  speed,  the  tool 
velocity  along  the  spiral  is  constant  at  v0  and  it  can  be  shown  that  the  time  t  is  given  by  the  path 
length  in  Eq.  (3)  divided  by  the  cutting  velocity  as  in  Eq.  (4), 

'(«)  =  d^C/S"  Eq.  (4) 

This  expression  for  the  relationship  between  0  and  t  in  Eq.  (4)  can  be  used  to  provide  the  axial 
position  of  the  tool  relative  to  the  workpiece  as  a  function  of  angular  position  0.  The  motion  of 
the  tool  in  the  axial  direction  (Z-direction)  is  given  by  Z=Asin(cot).  Given  the  effect  of  a  general 
offset  c,  as  described  in  the  preceding  section,  the  axial  position  is  modified  as  in  Eq.  (5), 

Z(t)  =  A(c  +  sin(2  *  n  *  fm  *  t))  Eq.  (5) 

Substitution  of  Eq.  (4)  in  Eq.  (5)  provides  the  position  of  the  tool  apex  in  polar  coordinates  (r,  0, 
z)  as  a  function  of  time  t  and  the  controllable  process  parameters  (fm,  A,  s0)  and  is  provided  in  Eq. 
(6), 


Z(0)  =  A(c  +  sin  (j&E.  j^R/s°  V02  +  I  de))  Eq.  (6) 

While  Eq.  (6)  describes  the  relative  tool  position  in  time,  the  actual  shape  left  by  the  tool 
interaction  with  the  work  surface  is  modeled  by  determining  the  relative  intersection  between  the 
work  surface  and  the  corresponding  tool  geometry  along  this  path.  In  this  case,  the  tool  geometry 
can  be  modeled  ideally  as  a  circle  or  sphere  of  radius  Rt  whose  instantaneous  position  is  defined 
by  Eq.  (6).  Boolean  subtraction  of  this  shape  from  the  work  surface  can  be  implemented 
numerically. 


A  similar  geometric  model  can  be  developed  for  sliding-type  texturing.  In  turning  configurations, 
the  motion  of  the  tool  takes  on  a  helical  path  across  the  outer  diameter  of  the  workpiece  at  a 
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constant  radial  position  (i.e.,  r  is  constant)  and  within  0-Z  space  as  in  Fig.  9(b).  At  constant 
spindle  frequency  (i.e.,  fw  is  constant),  the  angular  location  of  the  tool  is  defined  by  Eq.  (7), 

0(£)  =  2  *  7T  *  fw  *  t  Eq.  (7) 

Tool  motion  in  the  axial  direction  (i.e.,  Z-direction)  is  given  by  the  instantaneous  feed  of  the  tool 
and  the  superimposed  modulation  Asin(cot)  as  in  Eq.  (8);  substitution  of  Eq.  (7)  yields  an 
expression  in  terms  of  angular  location  in  Eq.  (9), 

Z(t)  =  s0  *  fw  *  t  +  A  *  sin  (2  *  n  *  fm  *  t)  Eq.  (8) 

z(0)  =  S^+A*s™(hf1)  £q-(9) 

Eq.  (9)  and  a  constant  value  for  the  radial  coordinate  describe  the  relative  tool  position.  The  shape 
of  the  work  surface  is  determined  by  solving  for  the  intersection  with  the  corresponding  tool 
geometry  along  this  path. 

5.2.2  Micro-Textured  Surface  Generation 

A  CNC  lathe  was  retrofitted  with  an  electro-mechanical  tool  holder  (M4  Sciences)  that  facilitated 
controlled  tool  motions  superimposed  on  the  machine  feed  in  the  axial  direction  of  the  workpiece, 
as  shown  in  Fig.  10(a).  A  polycrystalline  (PCD)  diamond  insert  (Hertel  HT-420CD)  with  a  nose 
radius  of  0.8  mm  was  used  for  finishing  the  work  surface.  A  carbide  grooving  insert  (Iscar  MIFR 
8-1.60-0.8)  with  a  nose  radius  of  0.8  mm  was  used  for  the  plunging-type  texturing  experiments 
and  a  carbide  threading  insert  (Iscar  MITR  8-MT  1-0.05)  with  a  nose  radius  of  0.05  mm  was  used 
for  sliding-type  texturing;  both  were  held  in  a  boring  bar  (Iscar  MG08-STFPR-09X)  attached  to 
the  electro-mechanical  tool  holder.  Tool  offsets  were  determined  by  white  light  interferometry 
(Zygo)  of  conventionally  grooved  workpiece  samples.  The  sinusoidal  waveform  (voltage, 
frequency)  to  the  piezoelectric  actuator  was  provided  using  an  electric  controller/amplifier  (M4 
Sciences).  Tool  displacement  was  monitored  using  a  capacitance  probe  (Capacitec).  Micro-scale 
surface  features  were  characterized  using  optical  microscopy  and  white  light  interferometry.  A 
Nikon  stereoscope  (SMZ800)  was  used  to  scan  the  surface  while  a  white  light  interferometer  was 
used  to  measure  the  three-dimensional  character  of  the  surface  features. 


(a) 


(b) 


Figure  10.  (a)  Experimental  machining  setup  and  (b)  ball-on-flat  wear  test  configuration. 


5.2.3  Wear  Testing  and  Coupled  Microstructure  Response 

Wear  experiments  were  carried  out  using  a  ball-on-flat  reciprocating  test  machine,  as  depicted  in 
Fig.  10(b).  Textured  and  non-textured  (control)  disks  were  tested  against  AISI  440C  stainless 
steel  balls  of  3  mm  radius  and  690  HV  hardness.  The  sliding  time  (20  min),  amplitude  (10.5  mm), 
frequency  (1.5  Hz)  and  normal  load  (34  N)  were  kept  constant  for  all  wear  tests.  2  mL  of  a 
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synthetic  poly-alpha-olefin  oil  (Synton  PAO  40)  were  applied  on  the  surfaces  prior  to  each  test 
and  no  additional  lubricant  was  added  during  the  tests.  Wear  performance  was  obtained  after 
three  tests  under  the  same  condition.  Volume  loss  (Vf)  was  determined  by  image  analysis  after  45 
wear  track  width  (Wt)  measurements  for  each  test.  Subsurface  microstructure  was  characterized 
using  electron  backscatter  diffraction  (EBSD)  with  a  Tescan  Mira3  XM  FEG-SEM  and  TSL  OIM 
analysis  software  to  obtain  crystallographic  orientation  information.  Samples  for  EBSD  were 
prepared  using  standard  mechanical  polishing  methods. 

The  effect  of  a  distribution  of  dimple  textures  on  wear  performance  was  measured.  Initial 
untextured  brass  and  aluminum  alloy  surfaces  were  produced  with  roughness  Ra  <  0.150  mm. 
Textured  surfaces  were  then  produced  using  these  workpieces  at  constant  spindle  speed  of  v0  = 
0.9  m/min,  modulation  frequency  of  fm  =  100  Hz  and  modulation  amplitude  of  A=0.015  mm. 
Figure  1 1  shows  optical  and  interferometric  images  of  the  plunging  type  textures  produced  in  the 
brass  and  aluminum  workpieces.  From  the  figure,  the  dimples  have  uniform  geometry  with  no 
discernable  differences  observed  between  the  work  materials.  Table  1  summarizes  the  dimple 
depth  measurements  obtained  from  interferometry.  From  the  table,  the  direct- write  texturing 
process  was  capable  of  generating  dimples  with  sub-micron  errors  in  depth  that  may  be  caused  by 
deformation  and/or  machine-tool  compliance.  Further,  the  surface  profilometry  in  Fig.  11  reveals 
a  negligible  region  of  material  pile-up  at  the  tool  exit. 


Table  1.  Experimental  surface  texturing  conditions  and  measurements  for  wear  tests. 
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Figure  11.  Surfaces  created  by  plunging-type  texturing  in  (a)  brass  360  and  (b)  AA6061-T6  at  A  =  0.02  mm, 
rt=  2.5  mm,  d  =  0  mm,  v0  =  12  m/min,  s0  =  2.0  mm/rev  and  fm  -  100  Hz. 

Tribological  wear  tests  were  carried  out  for  both  untextured  control  surfaces  (CS)  and  textured 
surfaces  (TS)  and  the  results  are  summarized  in  Fig.  12.  Under  the  experimental  conditions 
studied,  textured  surfaces  reduced  wear  with  respect  to  the  untextured  control  surfaces  on  both 
materials.  A  wear  reduction  of  32%  and  25%  was  achieved  by  texturing  the  surfaces  of  the  brass 
and  the  aluminium,  respectively.  Figure  12  shows  representative  wear  track  profiles  on  the 
textured  and  untextured  samples.  From  the  figure,  it  can  be  seen  that  the  wear  profiles  on  samples 
without  the  surface  textures  were  wider  and  deeper  than  those  on  textured  surfaces.  As  can  also 
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be  seen  in  Fig.  12(a),  the  amount  of  material  plastically  deformed  at  the  edges  of  the  wear  track 
was  also  reduced  for  the  textured  brass  samples. 
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Figure  12.  Wear  factor  measurements  (left)  for  control  surface  samples  and  textured  surface  samples  in 
brass  and  aluminum.  Wear  track  profiles  (right)  on  textured  (TS)  and  untextured  (CS)  samples. 


It  is  well  known  that  conventional  machining  can  induce  grain  refinement  and  phase 
transformations  in  the  work  subsurface  arising  from  severe  plastic  deformation  that  occurs  during 
the  chip  formation  process  [11,  12].  It  has  been  shown  elsewhere  [13]  that  machining  parameters 
of  undeformed  chip  thickness  and  rake  angle  have  primary  influence  on  spatial  extent  and 
severity  of  the  deformation  in  the  work  surface,  respectively.  In  plunging-type  texturing 
configurations,  these  parameters  vary  continuously  throughout  the  surface  generation  process.  To 
understand  the  effects  of  the  plunging-type  surface  texturing  process  on  surface  and  subsurface 
integrity,  a  controlled  experiment  was  conducted  to  characterize  the  microstructure  underneath  a 
model  dimple-type  surface  texture.  Sinusoidal  topographical  features  were  created  with  a  high¬ 
speed  steel  tool  with  a  rake  angle  of  a  =  0°,  the  modulation  amplitude  was  A  =  0.010  mm  and  the 
modulation  frequency  was  restricted  to  fm  =  1  Hz  to  facilitate  imaging  of  microstructure  using 
OIM  in  EBSD.  The  resulting  aspect  ratio  of  the  surface  textures  (i.e.,  2A/k)  was  0.2. 


The  microstructure  in  the  textured  subsurface  is  provided  in  Fig.  13.  From  the  figure,  the 
subsurface  micro  structure  featured  spatial  gradation  in  terms  of  grain  size,  width  of  the  severely 
sheared  zone  and  crystallographic  texture.  Broadly,  this  gradation  involved  a  saturated 
microstructure  field  featuring  an  equiaxed  ultra-fine  grained  (UFG)  morphology  on  the  exit  side 
of  the  sinusoidal  tool  motion  and  a  clearly  microcrystalline  structure  on  the  entry  side  of  the 
sinusoidal  motion.  The  extent  of  the  micro  structure  refinement  can  be  quantified  by  identifying 
grains  exhibiting  high-angle  misorientation,  generally  in  excess  of  15°.  The  average  grain 
diameter  where  grain  boundaries  exceed  this  threshold  was  measured  using  the  line  intercept 
method  and  is  given  by  615o  in  the  ensuing.  For  the  exit  side  of  the  sinusoidal  tool  motion  (left  in 
Fig.  13),  the  micro  structure  could  be  characterized  by  615o  =  0.5  pmover  a  thickness  of  5  pmin 
the  subsurface.  This  can  be  compared  to  the  entry  side  of  the  sinusoidal  tool  motion  (right  in  Fig. 
13)  where  615o  >  2  pm  is  observed  to  exist  over  the  subsurface  depth.  In  this  regard,  an 
asymmetric  micro  structure  field  comprising  of  UFG  domains  on  the  exit  side  and  coarse  grains 
on  the  entry  side  of  the  surface  texture  is  symptomatic  of  severe  shear  strains  being  imposed  in 
the  former  [14].  Crystallographic  textures  on  the  topographically  textured  surfaces  were  extracted 
from  these  OIM  graphs  and  revealed  a  similar  character  as  that  seen  in  conventional  orthogonal 
machining  with  constant  s0.  Texture  in  conventional  machined  surfaces  involves  concentration  of 
crystallographic  orientations  along  3  distinct  fibers  attributable  to  simple  shear  and  rolling-type 
deformation  geometries  [15].  However,  strengths  of  zones  within  these  fibers  were  seen  to  vary 
periodically  with  s0  during  topographic  texturing. 
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Figure  13.  Subsurface  micro  structure  and  texture  in  topographical  surface  texturing. 

The  gradation  of  microstructure  in  Fig.  13  is  a  consequence  of  the  periodically  varying  effective 
rake  angle  aeff(t)  that  manifests  in  variable  strains  and  depths  of  severely  sheared  zones  under  the 
generated  surface  [13].  The  effective  rake  angle  assumes  a  strongly  negative  value  in  the  exit  face 
of  the  tool  at  aeff  =  -32°  and  an  equivalent  positive  value  during  entry  of  the  tool  at  aeff  =  32°. 
From  previous  investigations  [15,  16],  it  is  clear  that  variation  in  a  from  -32°  to  32°  can  have 
profound  implications  on  deformation  and  corresponding  micro  structural  response  in  the  work 
subsurface. 

6.  Future  Studies  and  Implementation  on  Continuous  Processing  Systems 

The  above  results  in  development  of  integrated  severe  surface  plastic  deformation  (S2PD)  and 
modulation-assisted  texturing  (MAT)  platforms  provide  a  first-pass  attempt  to  link  microstructure 
and  texture  predictive  capability  in  FE  formulations  coupled  with  the  VPSC  framework.  From  the 
results,  wear  resistance  of  textured  surfaces  and  texture  evolution  within  surface  severely 
deformed  surfaces  can  be  obtained  by  a  calibrated  empirical  model  of  wear  given  an  initial 
surface  state  /  condition.  Calibration  of  the  FE  model  and  the  VPSC  parameters  was  enabled  by 
high-resolution  in  situ  digital  image  correlation,  as  well  as  empirical  measurements  of  the  as- 
deformed  micro  structure  (and  texture)  in  a  unit  model  impression  process.  Extension  of  the  unit 
process  in  the  FE  model  enabled  prediction  of  deformation  in  surface  mechanical  attrition. 
Coupled  effects  on  grain  refinement  introduced  by  the  texturing  process,  observed  using  electron 
microscopy  and  explained  due  to  spatially-varying  processing  parameters,  must  be  accounted  for 
in  the  generation  of  the  final  surface  by  the  hybrid  S2PD  and  MAT  processing  route. 

Toward  developing  a  fully-functional  surface  modeling  capability  to  determine  optimal  multi¬ 
scale  textured  surfaces  for  performance  applications,  the  integrated  modeling  capability 
developed  in  the  course  of  the  present  study  must  be  linked  with  calibrated  predictive  tools  for 
microstructure  variables  beyond  that  developed  here  for  crystallographic  texture  (e.g.,  grain  size, 
phase  transformations,  twinning)  wherein  the  deformation  conditions  in  both  the  micro  structure 
surface  generation  and  texturing  process  can  be  used  to  elucidate  the  final  surface  state.  A 
potential  approach  to  a  more  expansive  surface  modeling  capability  for  multi-scale  textured 
surfaces  could  be  found  in  leveraging  the  present  models/codes  with  extant  process- 
microstructure  (and  performance)  libraries  to  provide  first  pass  estimation  of  microstructure 
characteristics  based  on  the  wealth  of  information  regarding  grain  refinement  mechanisms  in  bulk 
SPD  studies.  Expansive  literature  on  grain  refinement  mechanisms  in  bulk  SPD  provide  a 
potentially  attractive  resource  in  this  regard  and  will  be  used  in  future  studies  to  establish  surface 
performance  models  of  significantly  higher  fidelity  and  overall  utility.  The  ultimate  goal  of  this 
area  of  research  is  to  facilitate  predictive  tools  to  a  priori  determine  processing  configurations 
needed  to  generate  plastically  graded  component  surfaces  with  custom-tailored  microstructures 
(e.g.,  twinned,  bimodal)  and  mechanical  properties. 
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